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(N 1 ABSTRACT 

We have made the largest-volume measurement to date of the transition to large-scale 
homogeneity in the distribution of galaxies. We use the WiggleZ survey, a spectroscopic 
survey of over 200 000 blue galaxies in a cosmic volume of ~ 1 /i~ 3 Gpc 3 . A new 
y\ . method of defining the 'homogeneity scale' is presented, which is more robust than 

methods previously used in the literature, and which can be easily compared between 
different surveys. Due to the large cosmic depth of WiggleZ (up to z = 1) we are 
able to make the first measurement of the transition to homogeneity over a range of 
cosmic epochs. The mean number of galaxies iV(< r) in spheres of comoving radius 
r is proportional to r 3 within 1 per cent, or equivalcntly the fractal dimension of 
the sample is within 1 per cent of D2 — 3, at radii larger than 71 ± 8 /i _1 Mpc at 
z ~ 0.2, 70 ± 5 ft- 1 Mpc at z - 0.4, 81 ± 5 h^Mpc at z - 0.6, and 75 ± 4 h^Mpc at 
z ~ 0.8. We demonstrate the robustness of our results against selection function effects, 
using a ACDM TV-body simulation and a suite of inhomogeneous fractal distributions. 
The results are in excellent agreement with both the ACDM W-body simulation and 
an analytical ACDM prediction. We can exclude a fractal distribution with fractal 
dimension below D2 = 2.97 on scales from ~ 80 /i _1 Mpc up to the largest scales 
probed by our measurement, ^300 /i~ 1 Mpc, at 99.99 per cent confidence. 

Key words: surveys - galaxies: statistics - cosmology: observations - large-scale 
structure of Universe. 
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1 INTRODUCTION 

One of the main assumptions of the standard theory of cos- 
mology, ACDM (based on cold dark matter and a cosmo- 
logical constant), is that the Universe is homogeneous and 
isotropic on large scales, and hence can be described by 
the Friedmann-Robertson- Walker (FRW) metric. 'Homoge- 
neous' means that its statistical properties (such as density) 
are translationally invariant; 'isotropic' means it should be 
rotationally invariant. The Universe clearly deviates from 
this on small scales, where galaxies are clustered, but on 
large enough scales (> 100/i - Mpc in ACDM), the distribu- 
tion of matter is assumed to be 'statistically homogeneous' 
- i.e., the small-scale inhomogeneities can be considered as 
perturbations, which have a statistical distribution that is 
independent of position. However, this is merely an assump- 
tion, and it is important for it to be accurately verified by 
observation. Over the last decade there has been a debate 
in the literature as to whether the Universe really is homo- 
geneous, or whether it has a fractal-like structure extending 
to large scales. It is important to resolve this contention if 
we are to be justified in assuming the FRW metric. 

In fact, although ACDM is based on the assumption 
of large-scale homogeneity and an FRW metric, inflation 
(which ACDM incorporates) actually predicts a certain level 
of density fluctuations on all scales. Inflation predicts that 
the primordial density power spectrum was close to scale- 
invariant. In the standard model, the scalar index n s , which 
quantifies the scale-depen dence of the pr imordial power 
spectrum, is close to 0.96 l|Baumannl 120091 ) . while a scale- 
invariant power spectrum has n 3 — 1. In this case, these 
density fluctuations induce fluctuations in the metric, 
which are virtually in dependent of s cale, and are on the or- 
der of S$/c 2 ~ 10" 5 (|Peacocklll999l ). Since these perturba- 
tions are small, the FRW metric is still valid, but it means 
that we expect the Universe to have a gradual approach to 
large-scale homogeneity rather than a sudden transition. 

The most important implication of inhomogeneity is the 
so-called 'averaging problem' in General Relativity (GR). 
This arises when we measure 'average' quantities (such as 
the correlation function and power spectrum, and parame- 
ters such as the Hubble constant) over a spatial volume. In 
doing so we assume the volume is homogeneous and smooth, 
when it may not be. Since the Einstein equations are nonlin- 
ear, density fluctuations can affect the evolution of the av- 
erage properties of the vol ume - th i s is kn o wn as the 'back- 
reacti o n mechanism' (e.g. Buchert J200C ) , Ellis 



e.g. Schwardl2002l; Kolb et alj|2005l; 


Rasanenll2006l; 1 Wiltshire! 


l2007al; BuchertJl200Sl: Rasanen 2011 


). This effect appears to 



Schward j2007f n" see iRasanenl (|201ll ) for a sum 



Buchertl 



|200d ). [Lr&: Scnwarzl l|2UUYI I; see 
mary). If we observe a quantity within such a volume, we 
need to take averaging into account to compare it with the- 
ory. It is therefore important to know how much inhomo- 
geneity is present, in order to obtain meaningful results from 
averaged measurements (and so most, if not all, cosmological 
measurements) . 

Backreaction has also been proposed as an explanation 
of dark energy, which is believed to be a negative-pressure 
component of the Universe that drives the accelerated ex- 
pansion. Some authors have suggested that instead of in- 
troducing exotic new forms of dark energy, or modifications 
to GR, we should revisit the fundamental assumptions of 
the ACDM model, such as homogeneity. If we assume that 
GR holds, but take inhomogeneities into account, it can be 
shown that backreaction can cause a global cosmic acceler- 



be too small to explain the observed acceleration, but high- 
lights the importance of understanding the amount of inho- 
mogeneity in the Universe. 

Another, related, consequence of inhomogeneity is that 
it can affect the path travelled by light rays, and the calibra- 
tion of clocks and rods of observers. It can therefore affect 
distanc e measurements, such as redshifts and luminosity dis- 
tances (|Wiltshirdl2009l ; lMeures fc Brunill2012l ). This has also 
been proposed as a possible explanation of the observed cos- 
mic acceleration, although the effec t appears to be on th e 
order of only a few percent at z ~ 1 (|Brouzakis et al]|2007l ) . 

Homogeneity is required by several important statisti- 
cal probes of cosmology, such as the galaxy power spectrum 
and n-point correlation functions, in order for them to be 
meaningful. Applying these to a galaxy sample below the 
scale of homogeneity would be problematic, since if a dis- 
tribution has no transition to homogeneity, it does not have 
a defined mean density, which is required to calculate and 
interpret these statistics. It is also not possible to model its 
cosmic variance, so the error in the measurements would be 
ill-defined, making it impossible to relate these statistics to 
a theoretical model. It is therefore important to quantify 
the scale on which the Universe becomes close enough to 
homogeneous to justify their use. 

Large-scale homogeneity is already well supported by 
a number of different observation s. In particular, th e high 
degree of isotropy of the CMB l|Fixsen et al.l Il996l ) gives 
very strong support for large-scale homogeneity in the 
early Universe, at redshift z ~ 1100. The isotropy of the 
CMB also indicates the Universe has remained homoge- 
neous, since there are no significant Integrated Sachs- Wolfe 
ISW) effects d istorting our view of the isotropic CMB 
Wu et al] Il999l ) . Other high-redshift evidence for homo- 
geneity includes the i sotropy of the X- Ray Background 
(XRB) ^Peebles! fl99^ ; IScharf et all l2000h . believed to be 
emitted by high- redshift sources, and the isotropy of radio 
sources at z ~ 1 (|Blake fc Wallll2002l ). 

However, these measurements of high-redshift isotropy 
do not necessarily imply homogeneity of the present Uni- 
verse. If every point in the Universe is isotropic, then this 
implies the Universe is homogeneous; so if we accept the 
Copernican principle that our location is non-spec ial, then 
the o bserved isotropy should imply homogeneity (|Peacockl 
ll999l )F1 However most of these measurements (except the 
ISW effect) only tell us about the high-redshift Universe. 
We know that it has evolved to a clustered distribution 
since then, and it is possible that it could also have be- 
come anisotropic. It is also possible for the matter distribu- 
tion to be homogeneous while the galaxy distribution is not, 
since the galaxy d istribution is biased relative to the matter 
field (H aiserlll984T ) - alth ough since galaxy bias is known to 
be lin ear on large scales (|Coleslll993l ; IScherrer fc ' Weinberg! 



1998 ), this seems unlikely. The ISW effect | |Sachs fc Wolfe 



19671 ) gives information about the low-redshift Universe, 



since it mostly probes the dark energy dominated era, z < 1 



1 We note that the Copernican principle is not incompatible with 
an inhomogeneous Universe. It assumes only t hat our location is 
non-special, not that every location is the sa me l ljovce et al .120001; 
IClifton et aT]|2008l ; ISvlos Labini et aljlioosh . 
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(jAfshordill200i ). but it is an integral over the line-of-sight 
and so does not give full 3D information. 

Galaxy surveys are the only 3D probe of homogene- 
ity in the nearby Universe, and a number of homogene- 
ity analyses have been carried out with different surveys, 
with seemingly conflicting results. Most statistical methods 
used to measure homogeneity have been based on the sim- 
ple 'counts-in-spheres' measurement, that is, the number of 
galaxies N(< r) in spheres of radius r centred on galaxies, 
averaged over a large number of such spheres. This quantity 
scales in proportion to the volume (r 3 ) for a homogeneous 
distribution, and homogenei ty is said to be re ached at the 
scale above which this holds. iHogg et al.l (|2005l ) applied this 
to the Sloan Digital Sky Survey (SDSS) luminous red galaxy 
(LRG) sample at z ~ 0.3, and found the data became con- 
sistent with homogeneity at ~ 70fe _1 Mpc for this galaxy 
population (using a different method of determining the ho- 
mogeneity scale than we do). 

This measurement can also be extended to a fractal 
analysis. Fractal dimensions can be used to quantify clus- 
tering; they quantify the scaling of different moments of 
galaxy counts in spheres, which in turn are related to the 
n-point correlation functions. The most commonly used 
is the correlation dimension D^ir), which quantifies the 
scaling of the 2-point correlation function, and is based 
on the counts-in-spheres, which scale as N(< r) ~ r° 2 . 
One can also consider the more general dimensions D q , 
where q are different moments of the counts-in-spheres. Us- 
ing fractal analyses, some researchers have found a tran- 
sition frcmi_D2_<C_3_to D% = 3, at around 70 - 150 ft. -1 



Mpc ( 


Martinez & Coles 


19941; Guzzo! 19971; Martinez et all 


19981; 


Scaramella et al.l 


19981; Amendola & Palladino! 1 19991; 


Pan & Coles! 2000 


; Kurokawa et al. 2001; Yadav et al.l 


20051; ISarkar et al. 


l2009h, whereas other authors have 



found no such transition (IColeman fc Pietronerol 1 19921; 



Pietronero et al.ll 19971: ISvlos Labini et a.1.1 119981: iJovce et al.l 
19991 ; ISvlos Labini et al.ll2009l ; ISvlos Labinill201lh . HoweverT 



many of the galaxy redshift surveys used in the above- 
mentioned works are too shallow, sparse, or have survey 
geometries too complicated, to give conclusive results. 

In this work, we use t he WiggleZ Dark Energy Sur- 
vey (|Drinkwater et aTll201C| l to make a new measurement 
of the counts-in-spheres and correlation dimension, to test 
for the transition to homogeneity. WiggleZ provides a larger 
volume than previous surveys, making it ideal for a homo- 
geneity measurement, and it covers a higher redshift range, 
allowing us to also investigate how homogeneity changes 
with cosmic epoch. It is not volume-limited, but we show 
that this does not significantly affect our measurement. The 
transition to homogeneity can be used as a test of a par- 
ticular cosmological model, since we would expect it to dif- 
fer for different cosmologies. In this work we test a ACDM 
model with best-fitting parameters from th e Wilkinson Mi- 
crow ave Anisotropy Probe (WMAP) data (|Komatsu et al.l 
120111 ). which we refer to as ACDM+WMAP. We demon- 
strate the robustness of our measurement against systematic 
effects of the survey geometry and selection function, by re- 
peating our analysis on both the GiggleZ A-body simulation 
and on a suite of inhomogeneous, fractal distributions. 

Before we make any meaningful test of homogeneity, 
however, it is crucial to properly define what we mean by 
the so-called 'scale' of homogeneity. Since there is only a 
gradual approach to homogeneity, such a definition may 



be arbitrary. In the past, authors have defined the 'scale 
of homogeneity' as the scale where the d ata becomes con- 
sistent with homogeneity within la (e.g. IHogg et a.1.1 120051 : 
iBagla et al]|200St lYadav et alj|2010l ). However, this method 
has several disadvantages. It depends on the size of the 
error bars on the data, and hence on the survey size. A 
larger survey should have smaller error bars, and so will au- 
tomatically measure a larger scale of homogeneity. It also 
depends on the bin spacing, and is susceptible to noise be- 
tween data points. We therefore introduce a different, and 
more robust, method for determining homogeneity: we fit a 
smooth, model-independent polynomial curve to all the data 
points, and find where this intercepts chosen values close to 
homogeneity. 

Certain parts of our analysis require the assumption 
of a cosmological model and, implicitly, homogeneity (i.e. 
for converting WiggleZ redshifts to distances, correcting 
for the selection function, calculating the uncertainties 
using lognormal realisations, and finding the best-fitting 
bias). In these cases, we use an input ACDM cosmology 
with h = 0.71, n m = 0.27, ST2 A = 0.73, D. b = 0.04482, 
erg = 0.8 and n B — 0.96. Here, the Hubble constant is 
Ho = lOO/i km s -1 Mpc -1 , fl m is the mass density, is 
the dark energy density, SI;, is the baryon density, erg is the 
root mean square mass variation within spheres of 8/i~ 1 Mpc 
radius, and n s is the spectral index of the primordial power 
spectrum. This is the same fiducial cosmology used by 
iBlake et al.l (|2011al ), and we use this for consistency. We 
discuss the implications of assuming a ACDM model on the 
results of our homogeneity measurement in Section [7] 

This paper is organised as follows. In Section [2] we de- 
scribe the WiggleZ survey and our dataset. In Section [3] 
we describe our methodology. We also explain our defini- 
tion of the 'scale of homogeneity' and present a new model- 
independent method for measuring this from data. In Sec- 
tion [4] we describe our analytic ACDM+WMAP model. We 
present our results in Section [5] We test the robustness of 
our method using fractal distributions and a ACDM A r -body 
simulation in Section [6] We discuss our results in Section [7] 
and conclude in Section [5] 



2 THE WIGGLEZ SURVEY 



The WiggleZ Dark Energy Survey jDrinkwater et al.|[201oh 
is a large-scale spectroscopic galaxy redshift survey con- 
ducted at the 3.9m Anglo- Australian Telescope, and was 
completed in January 2011. It maps a cosmic volume of ~ 1 
Gpc 3 up to redshift z = 1, and has obtained 239 000 redshifts 
for UV-selected emission-line galaxies with a median redshift 
of Zmcd = 0.6. Of these, 179 599 are in regions contiguous 
enough to be used for our analysis. It co vers ~ 1000 deg 2 
of equ atorial sky in 7 regions, shown in iDrinkwater et all 
(|2010T ) (their Fig.l). 

The observing strategy and gal axy selection criter ia of 
the WiggleZ survey are described in lBlake et al. | i|2009h and 
IDrinkwater et alj j201Cf) The selection function we use is 
described in IBlake et al.1 |2010l ). The targets were selected 
from Galaxy Evolution Explorer satellite (GALEX) obser- 
vations matched with ground-based optical measurements, 
and magnitude and colour cuts were applied to preferentially 
select blue, extremely luminous high-redshift star-forming 
galaxies with bright emission lines. 
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Table 2. Number of WigglcZ galaxies in each redshift slice. 



Redshift 


Number of galaxies 


0.1 < z < 0.3 


25187 


0.3 < z < 0.5 


45 698 


0.5 < z < 0.7 


70191 


0.7 < z < 0.9 


38 523 



The WiggleZ Survey offers several advantages for a new 
study of the scale of homogeneity. Its very large volume al- 
lows homogeneity to be probed on scales that have not pre- 
viously been possible, and at a higher redshift; it probes 
a volume at z > 0.5 comparable to the SDSS LRG cata- 
logue at z < 0.5. This allows us to make the first measure- 
ment of the change in the homogeneity scale over a range 
of cosmic epochs. We divide our sample into four redshift 
slices, 0.1 < 2 < 0.3, 0.3 < z < 0.5, 0.5 < z < 0.7 and 
0.7 < z < 0.9. The sizes of the WiggleZ regions in each red- 
shift slice are listed in Table[TJ they sample scales well above 
the expected scale of homogeneity. The numbers of galaxies 
in each redshift slice are listed in Table [2] 

We also benefit from having 7 regions distributed across 
the equatorial sky, which reduces the effect of cosmic vari- 
ance. In addition, WiggleZ probes blue galaxies, whereas 
SDSS (which has obtained the largest-scale measurements 
of homogeneity to date) used Luminous Red Galaxies, and 
so it can also constrain any systematic effects introduced 
by th e choice of tracer galaxy population (|Drinkwater et all 
I2010T ). Since blue galaxies are less biased, they are also more 
representative of the underlying matter distribution. 

There are, however, several aspects of the survey that 
could potentially be detrimental to a homogeneity measure- 
ment. WiggleZ has a complex window function, with a com- 
plicated edge geometry including holes in the angular cov- 
erage, and the spectroscopic completeness varies across the 
sky. In addition, the population properties of the galaxies are 
known to vary with redshift, due to the effects of Malmquist 
bias (since WiggleZ is a flux- limited survey), downsizing 
(the observed fact that the size of the most actively star- 
forming galaxies decre ases with time, ICowie et al.l I 19961 ; 
iGlazebrook et alj|2004l), and the colour and magnitude se- 
lection cuts (|Blake etal]|201Ch . This means that WiggleZ 
preferentially selects larger-mass, higher-luminosity galaxies 
at higher redshift. A consequence of this is that it is not 
possible to define volume-limited subsamples of WiggleZ. 
However, we can correct for these effects by using random 
catalogues (Se ction 13. II). which account for the survey selec- 



tion function (|Blake et al.ll20icih . We also divide the survey 
into four redshift-slices, reducing the amount of galaxy pop- 
ulation evolution in any region. In addition, we show that 
our results are not biased by the assumption of homogeneity 
in the selection function corrections, using an iV-body sim- 
ulation and a suite of inhomogeneous fractal distributions, 
described in Section [6] We are therefore confident that our 
result is not distorted by any features of the survey. 

We convert the redshifts of the WiggleZ galaxies to co- 
moving distances d c , using 



d c{z) = — 



where 



E(z) 



H(z) 
Ho 



dz' 
E(z'Y 



[ftm,0(l + z) 3 +OA,0] 1/2 , 



(1) 



(2) 



and we use the fiducial ACDM parameter values listed in 
Section 1. To do this, we assume the FRW metric and 
ACDM. This is necessary for any homogeneity measurement, 
since we must always assume a metric in order to interpret 
redshifts. Therefore in the strictest sense this can only be 
used as a consistency test of ACDM. However, if we find 
the trend towards homogeneity matches the trend predicted 
by ACDM, then this is a strong consistency check for the 
model and one that an inhomogeneous distribution would 
find difficult to mimic. We discuss this further in Section [7] 



3 METHODOLOGY 

Here we describe our methodology for measuring the transi- 
tion to homogeneity. We first calculate the mean counts-in- 
sphcrcs N(< r), then we find the fractal correlation dimen- 
sion D2(r), from the slope of N(r). Although they are closely 
related it is interesting to consider both, since the counts-in- 
spheres is the simplest measurement of homogeneity, whilst 
the correlation dimension provides direct information about 
the fractal properties of the distribution. We also describe 
our method of determining uncertainties, and our method 
of defining the 'homogeneity scale' Rh- 



3.1 Scaled counts-in-spheres A/"(< r) 

The simplest test of homogeneity of a set of points is to find 
the average number N(< r) of neighbouring points from any 
given point, up to a maximum distance r; if the distribution 
is homogeneous, then (for large enough r), 



N(< r) oc r 



(3) 



where D is the ambient dimension (the number of dimen- 
sions of the space; for a homogeneous volume, D = 3). 

We find N(< r) for spheres centred on each of the Wig- 
gleZ galaxies, and correct for incompleteness by dividing by 
the number expected for a homogeneous distribution with 
the same level of completeness. (We show that this does 
not bias our final results). This is done by finding the mean 
N(< r) about the coordinate position of the WiggleZ galaxy 
from 100 random catalogues, each with the same number 
density, window function and redshift distribution of the 
WiggleZ survey. The method of generatin g the random cat- 
alogues is described in lBlake et all (|2010r i. We then take the 
average over all the galaxies, to obtain the mean, scaled 
counts-in-spheres measurement Af(< r): 



(4) 



where G is the number of WiggleZ galaxies used as sphere 
centres, R is the number of random catalogues, N(< r) is the 
counts for WiggleZ galaxies, Nr(< r) the counts for random 
galaxies (centred on the position of the i th WiggleZ galaxy) , 
and pj = nw/n ra ,iid,i is the ratio of the total number of 
WiggleZ galaxies (nw) to the number of random galaxies in 
the j th random catalogue (ra ra nd.j)- In our analysis, we have 
G — nw, but this would not be the case if, for example, we 
excluded spheres near the survey edges. 

The random catalogue correction has the effect of re- 
ducing the scaling by the number of dimensions (i.e. D=3), 
so that for a homogeneous distribution Af(< r) scales as 
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Table 1. Comoving dimensions of WiggleZ regions in each redshift slice, in units of [h 1 Mpc], The dimensions correspond to the 
line-of-sight, RA, and Dec directions, respectively. 



Region 0.1 < z < 0.3 0.3 < z < 0.5 0.5 < z < 0.7 0.7 < z < 0.9 

00- hr 551 x 146 x 222 505 x 232 x 353 459 x 309 x 471 417 x 378 x 575 

01- hr 549 x 127 x 132 499 x 202 x 209 450 x 269 x 279 405 x 329 x 341 
03-hr 549 x 133 x 129 499 x 211 x 206 450 x 281 x 274 405 x 343 x 335 
09-hr 551 x 221 x 133 504 x 351 x 212 458 x 467 x 283 416 x 571 x 346 
11-hr 553 x 280 x 145 509 x 445 x 231 466 x 592 x 308 426 x 724 x 376 
15-hr 553 x 295 x 150 510 x 468 x 238 468 x 623 x 317 429 x 762 x 387 
22-hr 550 x 142 x 143 500 x 225 x 228 452 x 300 x 303 408 x 367 x 371 



10' 



10" 



10" 



WiggleZ 15-hr region, 0.5<z<0.7 



Mean volume of a shell surrounding 

a WiggleZ galaxy within the selection function 

True volume of shell 




r [h _1 Mpc] 



1000 



Figure 1. An illustration of the scales for which survey edge 
effects become important for the homogeneity measurement. Top 
panel: The mean volume V(r) in a thin shell of mean radius r 
surrounding a WiggleZ galaxy within the selection function (black 
curve). We show this for the 15-hr 0.5 < z < 0.7 region. The true 
volume of the shells is shown as a blue dashed line. The black 
curve deviates from the blue at large scales, where an increasing 
proportion of the shells extends outside the survey. Bottom panel: 
The ratio of the mean volume to the true volume, as a function 
of r. The grey dashed line indicates a ratio of 1. 



A/(< r) oc r 3 



(5) 



In each WiggleZ region we make A/"(< r) measurements 
in spheres with 12 to 15 logarithmically-spaced radial bins 
(depending on the size of the region). We determine the 
large-scale cutoff in each region by calculating the mean 
volume in the selection function, V(r), in a thin shell of 
mean radius r enclosing a central galaxy. We illustrate this 
for the 15-hr 0.5 < z < 0.7 region in Fig. fTJ and compare 
it with the 'true' volume of the shells. Below ~ 2/i~ 1 Mpc 
there is noise due to low resolution, but above this the two 
curves are very close. On large scales however, an increasing 
proportion of shells surrounding galaxies go off the edge of 
the survey, and their volume within the survey decreases. 
On such scales, corrections for edge effects will become im- 
portant. We take the large-scale cutoff of our homogeneity 
measurement at the radius of the maximum-volume shell. 
For the region in Fig. fT] this corresponds to shells that are 
~ 20 per cent complete. We show later that edge effects up 
to this scale do not impact our homogeneity measurement. 



Our correction method using random catalogues max- 
imises the use of the data, and accounts for incompleteness 
and the fact that WiggleZ is not volume-limited. However, 
it can potentially bias our result towards detecting homo- 
geneity, since it assumes homogeneity on the largest scales 
of the survey. It is equivalent to weighting each measure- 
ment by the volume of the sphere included within the sur- 
vey, multiplied by an arbitrary mean density. So while the 
counts-in-spheres measurement should have the advantage 
of not assuming a mean density (|Hogg et al.ll2005T ). our cor- 
rection method means that we do. Therefore A/(< r) should 
tend to 1 on the largest scales of the survey, regardless of 
whether homogeneity has been reached. However, for a dis- 
tribution with homogeneity size smaller than the survey, 
jV(< r) should reach 1, and remain at 1, for a range of scales 
smaller than the survey scale. We check the robustness of 
our method against effects of the selection function and cor- 
rection method in Section [6] and show that our analysis is 
robust against this potential source of systematic error out 
to scales far greater than the homogeneity scale we measure. 

There are other correction methods used in the litera- 
ture . ThereJs_tlie_s^c^Ued^ e^cd^ 

(e.g. IColeman fc Pietronerdll992l: |Pan fc Coles] |2000| . 120021 : 



ISvlos Labini et ai1l2009l ) which only considers central points 
that are surrounded by complete spheres within the survey. 
This therefore excludes as central points any galaxies within 
a certain distance from the survey edges. However, this does 
not make the best use of the data, since it excludes data and 
so reduces the volume of the sample. 

There is also the so-called 'angular correction' model, 
which has been shown to be more optimal, by using all the 
availa ble data without int roducing a bias due to edge correc- 
tions (|Pan fc Coles! 120021 ). This corrects measurements in a 
sphere by the solid angle subtended by regions in the sphere 
that are outside the survey boundary. However, as they point 
out, this method is difficult to apply to surveys with a com- 
plicated geometry, especially if they contain holes, as Wig- 
gleZ does. Finally, another correction that minimises bias at 
the surv ey edges, but wastes little data, is the 'Ripley' es- 
timator (|kiplevlll977l ;l Martmez et al.|[l99ct ). which corrects 
the measurement in a sphere by the area of the sphere con- 
tained within the survey. Both the angular correction and 
the Ripley estimator assume isotropy of the samples. Due 
to the geometry of the WiggleZ survey, we choose to use a 
random catalogue correction, but make robustness tests to 
quantify any bias it may introduce to the results. 

To demonstrate the robustness of our measurement 
against the method of correcting for the selection function, 
we compare our method to an analysis using only complete 
spheres, with and without correcting for incompleteness, in 
Section \6. 31 and show we obtain consistent results. 
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3.2 Correlation dimension Z?2(r) 

Fractal dimensions can be used to describe the clustering 
of a point distribution. There exists a general family of di- 
mensions D q , the Minkowski-Bouligand dimensions, which 
describe the scaling o f counts in spheres cen tred on points 
(see e.g. iBorgamll 19951 ; iMartmez fc Saarll2002l . for a review) . 
To completely characterise the clustering of our galaxy dis- 
tribution we would need to consider all the moments q of 
the distribution (corresponding to combinations of n-point 
correlation functions). However, to identify the scale of ho- 
mogeneity we consider only D2, the 'correlation dimension,' 
which quantifies the scaling behaviour of the two-point cor- 
relation function £(r). 

If we take a galaxy and count the number of other galax- 
ies, N(< r), within a distance r, then this quantity scales as 



N(< r) <x r 2 , 

where D2 is the fractal dimension of the distribution. 
From this, the correlation dimension is defined as: 



D 2 (r 



d\nN(< r) 
rflnr 



(6) 



(7) 



Since we must correct each WiggleZ iV(< r) mea- 
surement for completeness, obtaining the scaled quantity 
JV(< r) oc r D2 ~ 3 , we must calculate Dz(r) via 



omr 



(8) 



For a homogeneous distribution, D2 = 3. If D2 < 3 
then the distribution has a scale-invariant, fractal (and so 
inhomogeneous) clustering pattern. If D2 > 3 the distribu- 
tion is said to be 'super- homogeneous' and co rresponds to a 
lattice-like distribution (|Gabrielli et al.|[2002l ). A power-law 
in 1 + £(r) ~ r" 7 has D 2 = 3 - 7 for £ > 1. 

Some previous works have found Dzir) simply by fit- 
ting a straight line to a log- log plot of N{< r). This method 
can give a false indication of a fractal (M artinez et al.l 
Il998l ); calculating D2(r) explicitly gives a more reliable 
measurement. 



We note that our estimator for Af(< r), Eq.[H is essen- 
tially equivalent to 1 + f ff (r), where (|Hamilton| [l992) 



£(r) = 



x £(x)dx. 



(This can be seen more clearly by rearrangement of the 
theoretical expression for A/"(< r) given by Eq. I19|) . Many 
measurements of £ (r) have been made with different 



galaxy surveys (e.g . Hawkins et alj 20031 ; IZehavi et al.ll2005l : 



iBlake et al1l2011bl ; lBeutler et al.ll201ll ). On small scales the 
correlation function is well described by a power law, 



(10) 



where ro ~ 5/i _1 Mpc is the so-called clustering length, and 
7 w 1.8, depending on the galaxy population. On scales 
> 20h~ Mpc however, the correlation function is observed 
to turn over, consistent with large-scale homogeneity. 

However, the correlation function cannot be used to test 
large-scale homogeneity, since the way it is determined from 
surveys depends on the mean galaxy density . Determinations 
of £(r) commonly use the Landy-Szalay l|Landv fc Szalavl 



Il993h or Hamilton estimators l|Hamiltonlll993l ). which com- 
pare the galaxy clustering to that of random catalogues of 
the same mean density as the survey. 

Our A/"(< r) estimator also compares the data to ran- 
dom catalogues, since we must correct for the selection func- 
tion. Therefore J\f(< r), like £(r), does assume a mean den- 
sity on the scale of the survey. However, our estimator is 
slightly different, as we correct each object separately rather 
than an average pair count. 

The correlation dimension 1)2(7*) > on the- other hand, 
measures the scaling of J\f(< r), which is not affected by the 
assumption of the mean density. (This only affects the am- 
plitude of J\f(< r)). Deviations from a volume-limited sam- 
ple, which require random-catalogue corrections, only cause 
second-order changes to D2(r), whilst they would be leading- 
order in the raw correlation function. Therefore, Z?2(r) is 
much more robust to both the assumed mean density and 
details of the selection function, making it the most reliable 
measure of homogeneity. 



3.3 Lognormal realisations and covariance matrix 

Lognormal realisations (|Coles fc Jonesiri99ll ) are an impor- 
tant tool for determining uncertainties in galaxy surveys. A 
lognormal random field is a type of non-Gaussian random 
field, which can be used to model the statistical properties 
of the galaxy distribution, and simulate datasets with an 
input power spectrum. We have used 100 such realisations, 
generated using an input ACDM power spectrum with the 
fiducial parameters listed in Section [T] These are sampled 
with the survey selection function, to create 100 mock cat- 
alogues for each of the WiggleZ regions. We use these to 
calculate the full covariance and errors of our measurement. 
Jack-knife resampling does not permit enough independent 
regions within the survey volume to give a reliable estimate 
of the uncertainties. 

We obtain Af(< r) and D2(r) for each of the lognormal 
realisations in the same way as for the WiggleZ data. The 
covariance matrix between radial bins i and j is given by 



I " 

— $>-(n) - xMnxfa) - x(rj)], (11) 



/g\ reausaric 



where x(r) is N{< r) or D%(r), the sum is over lognormal 
realisations I, n is the total number of lognormal realisations 
and x(r) = — X^ILi x i( r )- The diagonal values j = k give the 
variance, a 2 . 

The correlation coefficient between bins i and j is given 

by 



yj CaCj 



(12) 



The correlation matrices for A/"(< r) and -D2M are 
shown in Figures [2] and [3] respectively, for the combined re- 
gions (see next section) in the 0.5 < z < 0.7 redshift slice. 

It is noticeable that the Af(< r) measurement is more 
correlated than Da(r). This is because A/"(< r) is effectively 
an integral of D^fr), so its covariance is effectively a cu- 
mulative sum of that of D2(r). It can also be explained by 
the fact that D2(r) is the logarithmic slope of Af(< r), so 
it does not depend on the correlations between widely sepa- 
rated A/"(< r) bins but rather the variations between neigh- 
bouring bins. 
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We note that the uncertainties calculated using lognor- 
mal realisations assume ACDM, and represent the variance 
we would expect to measure in a ACDM universe. It would 
not be possible to calculate uncertainties for a fractal uni- 
verse, since a fractal has no defined cosmic variance as it has 
no defined mean density. However, calculating uncertainties 
this way gives a valid consistency check of ACDM. 

3.4 Combining WiggleZ regions 

In each of the 4 redshift slices, we combine the measurements 
from each of the 7 WiggleZ regions using inverse- variance 
weighting. The combined measurements i mm b are given by 



■Ecomb I ^ J 



(13) 



where n are the 7 WiggleZ regions, x n (i) is the measurement 
in the i th radial bin in the n th region, and u>„(r) = 1/C n (i, i) 
is a weight function. C n (i,j) is the covariance matrix of the 
n region, and the combined covariance matrix is calculated 
by 

Ccomb(^; j) = 



En Cn(i, j)w n (j)w n (j) 
E„™4*)E„ W n(j) 



(14) 



We show C CO mb(i,i) for the A/"(< r) and D2(r) mea- 
surements in the 0.5 < z < 0.7 redshift slice, in terms of the 
correlation matrix dj / y/CuCjj , in Figs. [2] and [3] 



3.5 Likelihood (x 2 ) & parameter fitting 

The x 2 °f our model fit to the data is given by 



X 2 = [ 2; th(r-i)-2;ob s (r0]C 4 /[2;th(rj)-2; obs (r j )], (15) 

i=l j=i 

where x(r) is jV(< r) or D2(r), ribins is the number of radial 
bins, Xth is the theoretical model and :r b s is the measured 
value. CJ" 1 is the inverse of the covariance matrix. 



3.6 A new method of denning the 'homogeneity 
scale' Rh 

As previously mentioned, the definition of a 'homogene- 
ity scale' is somewhat arbitrary since we expect the Uni- 
verse to hav e a gra dual approach to h omog eneity. Recently 
iBagla et al] (120081 ) and lYadav et all l|2010f) proposed that 
the homogeneity scale be defined where the measured frac- 
tal dimension D q becomes consistent with the ambient di- 
mension within the lcr statistical uncertainty, UAD q - They 
derived an approximation for D q (r) and uad, , in the limit of 
weak clustering, for a given correlation function, and showed 
that both scale the same way, and so the homogeneity scale 
stays constant, with bias and epoch. This definition is there- 
fore beneficial as it is not arbitrary, and is robust to the 
tracer galaxy population. However, in deriving <tad„ they 
considered only shot noise, and cosmic variance from vari- 
ance in the correlation function, while ignoring contributions 
from the survey geometry and selection function. Any real 
survey will have these contributions to the statistical uncer- 
tainty, which cannot be separated from the variance due to 
the correlation function alone. The value derived from this 
definition in a real survey is therefore difficult to interpret 
in a meaningful way, i.e. one that allows comparison with 
theory, or with different surveys of differing volume and se- 
lection function. 

We therefore introduce a different method for deter- 
mining a 'homogeneity scale' Rh, which is easier to com- 
pare with theory and between surveys. Our method is to fit 
a smooth, model-independent polynomial to the data, and 
find the scale at which this intercepts a chosen value, or 
'threshold,' close to homogeneity. This scale is then defined 
as the homogeneity scale Rh - For example, Rh could be the 
value of r at which the polynomial intercepts a line 1 per 
cent from Af(< r) — 1, or D^ir) = 3 (see Fig. |4]for an illus- 
tration). The uncertainty is found using the 100 lognormal 
realisations. The homogeneity scale measured this way does 
not depend directly on the survey errors (although the un- 
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Figure 4. Illustration of our method of defining the homo- 
geneity scale, Rh, shown here for the D2(r) measurement. We 
first fit a model-independent polynomial (red curve) to the data 
(black data points). We then find where this intercepts a chosen 
value close to homogenenity, e.g. 1 per cent from homogeneity, 
D 2 = 2.97 (dotted grey line). This gives us Rh- We find the 
uncertainty in Rh from the root mean square variance of 100 
lognormal realisations (pink curves). 



certainties on Rh do), and is less susceptible to noise in the 
data, making this a preferable method that allows compar- 
isons between different surveys. It also allows easy compar- 
ison between the data and a given model, e.g. ACDM, and 
we can check that the data converges to N — 1 or D2 = 3 
as expected for a homogeneous distribution, by choosing a 
range of thresholds approaching homogeneity (see Fig. [5]). 

We can also take this further and construct a likelihood 
distribution for the homogeneity scale, as described in Sec- 
tion [53] 

Although our choice of D2 threshold is arbitrary, by 
choosing the same threshold for different surveys we obtain 
an Rh value that can be meaningfully compared, and that 
can be easily compared to a theoretical model. Our choice 
of threshold may be limited by the amount of noise in the 
data, however. For instance, we can measure an intercept 1 
per cent away from homogeneity for the WiggleZ data, but 
cannot measure 0.1 per cent in two of the redshift slices, 
due to noise (the data do not come this close to homogeneity, 
although they are consistent with it within the uncertainties) 
- see Fig. [9] Nonetheless, we can easily choose a threshold 
that is possible given our data, and use this to compare with 
a model. The Baryon Acoustic Oscillation (BAO) feature 
can also potentially affect the appropriate choice of intercept 
value as it causes a small distortion in D2(r) - we discuss 
this in Section [7] 

For the main results in this paper we have chosen a 
threshold of 1 per cent away from homogeneity, since this is 
about the closest threshold to homogeneity we can measure, 
considering the noisiness of the data. 



4 ACDM MODEL PREDICTION OF M{< R) 

AND D 2 (R) 

In this section we derive theoretical ACDM predictions for 
the counts-in-spheres and correlation dimension. This allows 
us to compare our measurements of the transition to ho- 



mogeneity to the predictions of a ACDM model that fits 
WMAP data. 



4.1 Af(< r) and D 2 (r) 

For a particular galaxy population, we can calculate the 
mean counts-in-spheres N(< r) from the 2-point matter cor- 
relation function predicted by ACDM. The 2-point correla- 
tion function £(r) is defined as the excess probability above 
random of finding tw o objects in vo lumes dVi and dVz, sep- 
arated by distance r (|PeebresHl98(ih : 



P(r) = p 2 [l + £(r)]dlW 2 . 



(16) 



where p is the mean number density. 

The mean number of galaxies surrounding a random 
galaxy up to distance r is found by integrating the correla- 
tion function 

N{< r)=p [ [1 + b 2 i{r')]A-Kr' 2 dr' , (17) 
Jo 

where b is the galaxy bias, relating the clustering of a partic- 
ular galaxy population to the underlying dark matter distri- 
bution. Note that ACDM assumes large-scale homogeneity, 
and indeed we must assume large-scale homogeneity in order 
for a mean density p to be defined. 

We obtain our model correlation function by transform- 
ing a ACDM matter power spe ctrum Pss(k) generated us- 
ing C AMB ijLewis et al.l l2000h . Since we make our mea- 
surements in redshift space, we first convert Pss(k) to the 
redshift-space galaxy power spectrum Pg{k) (described in 
the next section), then convert this to the redshift-space 
galaxy correlation function £ 9 (s), where s denotes distance 
in redshift-space. Since we use the angle-averaged power 
spectrum (assuming the power spectrum is isotropic), we 
do not need to integrate the angular part of the fe-space in- 
tegral, and so use a spherical Hankel transform rather than 
a Fourier transform to obtain £ 9 (s): 



^w = 2^/™^r* 3d *- (18) 

To compare with our WiggleZ measurement (Eq. Q, 
where we correct for incompleteness, we divide our counts- 
in-spheres prediction by the number that would be expected 
for a random distribution, i.e. p^irr 3 : 

3 



Af(< 



47rr 3 



[l + ^(s)]47rs / ds. 



(19) 



We calculate the model D%(r) by simply applying Eq. 
[7]to our model Af(< r). 

4.2 Redshift-space distortions and nonlinear 
velocity damping 

Here we describe how we implement redshift-space distor- 
tions in our analytical model. In practice, we measure the 
positions of galaxies in redshift-space, which are affected 
by redshift-space distortions. These are due to the peculiar 
velocities of galaxies along the line of sight, which add to 
the measured redshifts and perturb the inferred galaxy po- 
sitions. This anisotropic effect creates anisotropy in the ob- 
served redshift-space galaxy power spectrum Pg(k,p), and 
can be modelled by multiplying (convolving in configuration 
space) the real-space matter power spectrum by an angle- 
dependent function F(k,p): 
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Table 3. Values of the redshift space distortion parameter /3, 
and the pairwise velocity dispersion <r p , used in our modelling 
of nonlinear redshift-spa c e disto rtion effects. These values were 
obtained by iBIake et al,l l|2011ah from fits to the WiggleZ two- 
dimensional power spectrum. 



Redshift 





a p [h km/s] 


0.1 < z < 0.3 


0.69 


346 


0.3 < z < 0.5 


0.73 


275 


0.5 < z < 0.7 


0.60 


275 


0.7 < z < 0.9 


0.51 


86 



P 9 *(fc,ju) = b 2 F(k lf i)P 66 (k). 



(20) 



There are two forms of redshift-space distortion of rel- 
evance to our measurement, which we find are necessary for 
a good fit to the data: 

(i) On large, linear scales (> 20h~ Mpc) the bulk infall 
of galaxies towards overdensities creates an enhancement in 
the observed power spectrum along the line of sight. Thi s 
can be modelled by the linear Kaiser formula (|KaiserJll987h , 



P 9 s (fc, M ) 



}JL 2 
b 



(21) 



where Pss,h{k) is the linear matter power spectrum, / is the 
growth rate of structure and j3 = / jb is the redshift-space 
distortion parameter. 

Note that this formula assumes a perturbed FRW universe 
with small real-space density perturbations |<5(r)| <g 1. 

(ii) On quasilinear scales (10 < s < 20ft" 1 Mpc ) the pe- 
culiar velocities resulting from the scale-dependent growth 
of structure distort the shape of the power spectrum, via 
a scale-dependent damping effect. A common way o f mod- 
elling this is the 'streamin g model' l|Peeblesl Il980l ; iFisherl 
119951 ; lHatton fc Cold[l998h , which combines the linear the- 
ory Kaiser formula with a velocity streaming term. We 
choose to use a Lorentzian term, F = [1 + (kap/i) 2 ]^ 1 , for an 
exponential velocity probability distribution function, since 
IBIake et al.l (|2011al ) found this to be a good fit to the Wig- 
gleZ power spectrum for k < 0.1/t~ 1 Mpc. This give s us the 
so-called 'dispersion model' l|Peacock fc Doddslll994r i for the 
full redshift-space power spectrum: 



P s a {k^) = b' 



(1+/V) 2 

1 + (k(7pfl)' 



■Pss(k). 



(22) 



Here, a p is the pairwise velocity dispersion along the line 
of sight. Both dp and ft axe parameters that must be fit- 
ted to the data. We use the values obtained bv lBlake et al.l 
l)2011al ) in each redshift slice from fits to the WiggleZ two- 
dimensional power spectrum - these are listed in Table [3] 

We note that the streaming model is motivated by viri- 
alised motions of particles within haloes, on much smaller 
scales - the so-called 'Finger of God' effect at < 2h~ l Mpc. 
However, it is heuristic in natu re and can al s o descr ibe phys- 
ical scales of tens of /i -1 Mpc. IBIake et al.l l|2011aT ) apply it 
this way by fitting for a p on these scales, rather than on 
Finger-of-God scales. We find that including it gives a sig- 
nificant improvement of our model fit to data. 

To obtain the angle-averaged redshift-space galaxy 
power spectrum Pg(k) we need to convert the full F(k,n) 
to F(k), which we do by integrating over \i: 



F(k) 



1 (1 + /?M 2 ) 2 
=0 (1 + (fca. uM )2) 



dfj,. 



(23) 



The angle- averaged redshift-space galaxy power spec- 
trum is then 



P°(k) = b 2 F(k)P ss (k). 



(24) 



Eq. fJH normally assumes a linear matter power spec- 
trum; however, we choose to use a nonlinear Pss(k), cal- 
culat ed by CAMB using the H ALOFIT code |Smith et all 
2003), since IBIake et al" l|2011al ) found this gave a better fit 
to the data. 



4.3 Correcting the WiggleZ data for galaxy bias 

The amplitude of the galaxy correlation function is affected 
by galaxy bias and redshift-space distortions, and the shape 
is affected by nonlinear damping. Therefore, these also affect 
the amplitude and shape of A/"(< r), as well as the measured 
scale of homogeneity. It is possible to correct the data for 
bias, and so determine the homogeneity scale for the under- 
lying matter distribution, by assuming a particular model, 
i.e. our ACDM + WMAP model, fitting for bias by minimis- 
ing the x 2 value, and correcting the data for this. 

In our analysis, we only consider linear galaxy bias. This 
relates the galaxy correlation function £ 9 (r) to the matter 
correlation function £ m (r) through £ 9 (r) = 6 2 £ m (r). Since 
we are only interested in large scales, we do not consider 
scale-dependent bias that may occur on small scales. 

We fix the redshift space distortion parameter fi and 
velocity dispersion o p to the values listed in Table [3] We 
then obtain our corrected measurements A/biasfrec(< r), 

Af(< r) 



A"bia 



.(<r) = 



b* 



+ 1, 



(25) 



and we calculate the bias-corrected correlation dimension 
Ai.biasfrccM f rom tnis > using Eq. [7] 

Since we assume a ACDM model to fit for bias, we 
cannot make a model-independent measurement of transi- 
tion to homogeneity of the underlying matter distribution. 
However, it is still interesting to look at the variation of the 
homogeneity scale of the matter distribution with redshift, 
assuming ACDM + WMAP. 

We note that our measurement of the homogeneity 
scale of the WiggleZ galaxies is independent of the ACDM 
modelling shown in this section. This modelling is only done 
so that we can show that the measurement is consistent 
with that expected from ACDM. 



5 RESULTS 

5.1 Af(< r) and D 2 (r) 

The J\f(< r) measurements in each of the four redshift 
slices are shown in Fig. [5] The data are compared with a 
ACDM+WMAP model (described in Section gj). For each 
successive redshift slice, the reduced \ 2 values are 0.57, 0.91, 
0.69 and 1.1. The first two redshift slices have 14 data bins 
from 12.5 to 251 /i -1 Mpc, while the last two have 15 data 
bins from 12.5 to 316 ft _1 Mpc. The data is consistent with 
a monotonically decreasing function, so we can fit a polyno- 
mial and find where this intercepts a chosen threshold, as 
per our definition of homogeneity. 
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The intercepts of the polynomial fit with N = f .01 (1 
per cent away from homogeneity), which we define as the 
homogeneity scale Rh, are shown as red error bars. The 
errors were determined from the 100 lognormal realisations, 
and correspond to the square-roots of the diagonal elements 
of their covariance matrix. The Rh values and their errors 
are shown in Fig. [5] and listed in Table 2] along with the 
values for the ACDM model, which are in good agreement. 

The D2 (r) measurements in each of the redshift slices 
are shown in Fig. [5J along with a ACDM+WMAP model 
with best-fitting bias. In each redshift slice the data range 
and degrees of freedom are the same as for the Af(< r) mea- 
surement. The reduced x 2 values in each redshift slice are 
0.83, 0.90, 0.74 and 0.98. In each case a polynomial is fitted 
to the data. The homogeneity scales measured where these 
intercept 1 per cent away from homogeneity, D2 = 2.97, are 
listed in Table [4] and are also in excellent agreement with 
the ACDM values. 

5.2 Effect of bias and ag(z) on Rh 

The homogeneity scale of the model galaxy distribution, 
measured at 1 per cent from homogeneity, will depend on 
the amplitude and shape of the correlation function, and so 
on galaxy bias b and ag(z). We would also expect it to de- 
pend on redshift, since ag(z) increases over time. These two 
parameters are in fact completely degenerate in the A/"(< r) 
and D^(r) measurements. 

Fig. [3 shows how our ACDM Af(< r) and D 2 (r) models 
vary with bias, at z = 0.2, for fixed erg (2 = 0) = 0.8. Larger 
bias means a larger amplitude of clustering, so that both 
curves are steeper on small scales. This means that the mod- 
els reach homogeneity at larger radii for higher bias. This 
can be understood qualitatively, since highly biased galax- 
ies are more clustered together than less biased galaxies, so 
we must go to larger scales before we reach a homogeneous 
distribution. 

Fig. [5] shows how the homogeneity scale Rh of the 
galaxy distribution varies with bias, for different intercept 
values approaching homogeneity, for the A/"(< r) and D2(r) 
models at z = 0.2 with fixed cr$(z — 0) = 0.8. For a par- 
ticular intercept value, larger bias again gives a larger ho- 
mogeneity scale. Since Af(< r) and -D2M approach homo- 
geneity asymptotically, the intercept jumps to higher values 
when we consider an intercept value closer to homogeneity. 
This image also shows the mapping between A/"(< r) and 
D2(r) homogeneity values for intercepts at 1, 0.1 and 0.001 
percent away from homogeneity. They are not identical since 
they are slightly different methods, but give similar results. 
This plot illustrates that there are many potential ways to 
define homogeneity, and so it is important to make consis- 
tent measurements between surveys in order for them to be 
comparable with each other and with theory. 

Our bias-corrected Af(< r) and D2(r) measurements 
are listed in Table [4] The errors on the bias-corrected data 
were determined by applying a bias correction to each of 
the 100 lognormal realisations individually, and recalculat- 
ing the covariance matrix. Since the bias correction aims to 
set the bias of all the realisations to b = 1, it lowers the 
overall variance, and so the error bars are slightly smaller 
than for the pure data. These measurements give our mea- 
sured homogeneity scale for the matter distribution, assum- 
ing ACDM. We find that this scale increases with redshift, as 
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Figure 8. Rh for ACDM N{r) (green) and -D2M (purple) mod- 
els at z = 0.2 with differing bias b 2 . Each curve corresponds to 
Rh evaluated at a different threshold - 1, 0.1 or 0.01 per cent 
away from homogeneity (from bottom to top, labelled). 



expected in ACDM. However, since we are assuming ACDM, 
which has ag(z) increasing with time, this is not a model- 
independent result. 

As already mentioned, the effect of bias on the ho- 
mogeneity scale of galaxies is degenerate with the ampli- 
tude of the correlation function, <Ts(z), since the correlation 
function at redshift z depends on a combination of these, 
£(r, z) oc b 2 ag(z) 2 . So far we have assumed a fixed value of 
ag (z = 0). But we can also make predictions independent of 
ag, by finding how Rh changes as a function of the combina- 
tion b 2 ag(z) 2 . This is shown in Fig.||] We also show the Wig- 
gleZ results, which we have plotted for the best-fit b 2 ag(z) 2 
value in each redshift slice. This increases with redshift, so 
the data points from left to right go from low to high red- 
shift. The WiggleZ results are in very good agreement with 
the ACDM+WMAP predictions. We show, for comparison, 
the values obtained from defining the homogeneity scale as 
where the data comes within la of homogeneity. These val- 
ues have much greater stochasticity than those from our 
method of fitting a smooth curve to many data points, and 
do not give informative results in this plane. 

We see that the model RH-b 2 ag(z) 2 curves are mono- 
tonically increasing. Since we expect ag(z) in ACDM to grow 
over time due to growth of structure, we would therefore 
also expect the homogeneity scale to increase over time, for 
galaxies with fixed bias. 

For the WiggleZ data, however, the measured homo- 
geneity scale does not appear to decrease with redshift. This 
is explained by the fact that the WiggleZ galaxies have in- 
creasing bias with redshift, assuming a ACDM growth rate. 
As explained previously, this is understood to be due to the 
effects of Malmquist bias and downsizing on the selection of 
the WiggleZ galaxy population, and the colour and magni- 
tude cuts. This counteracts the effect of decreasing ag with 
redshift. As can be seen in Fig. [9] within the b 2 ag(z) 2 range 
of the data we would not expect a significant change in Rh , 
measured at 1 per cent from homogeneity, with redshift. 
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Figure 5. Scaled counts-in-sphcres A^(< r) for the combined WiggleZ data in each of the four redshift slices (black error bars). A ACDM 
model with best-fitting bias b 2 is shown in blue. A 5 th -degrec polynomial fit to the data is shown in red. The red errorbar and label 
show the homogeneity scale Rh for the galaxy distribution, measured by the intercept of the polynomial fit with 1.01 (1 per cent away 
from homogeneity), with the error given by lognormal realisations. This scale is consistent with the ACDM intercept with 1.01, labelled 
in blue. 



Table 4. Measured values of the homogeneity scale Rjj (where data intercepts 1 per cent of the homogeneity value). The Rfj values 
shown arc all for the galaxy distribution, except the values for the bias-corrected data, which arc for the underlying matter distribution, 
i.e. b 2 = 1. The bias-corrected values directly assume a ACDM+WMAP model. 
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5.3 Likelihood analysis for homogeneity scale 

Rather than trying to measure the scale of homogeneity di- 
rectly, it can be informative to consider the likelihood that 
the Universe has reached homogeneity by a certain scale. We 
can construct a probability distribution function (PDF) for 
the homogeneity scale, by combining the likelihood distri- 



bution of the data with that of the definition of the homo- 
geneity scale. This gives the probability that homogeneity is 
reached at a certain scale r. 

We could apply a likelihood analysis to either our 
J\f(< r) or -D2M measurements. However, it would arguably 
be invalid for highly correlated data, such as the 7V(< r) 
measurement, since the different contributions to the prob- 
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Figure 6. Same as for Fig. [5] but for the correlation dimension D^ir). The Diir) measurements for the combined WiggleZ data in 
each of the four redshift slices are shown as black error bars. A ACDM model with best-fitting bias b 2 is shown in blue. A 5th-dcgrcc 
polynomial fit to the data is shown in red. The red errorbar and label show the homogeneity scale Rh measured by the intercept of the 
polynomial fit with 2.97 (1 per cent away from homogeneity), with the error given by lognormal realisations. This scale is consistent with 
the ACDM intercept with 2.97, labelled in blue. 




Figure 7. The effect of bias on a ACDM A^(< r) model (left) and £)a(r) model (right) at z = 0.2. Increasing bias increases the value 
of J\f(< r) on small scales, and decreases the value of D2(r) on small scales, and produces a larger homogeneity scale, as seen by the 
intercepts of the curves with 1 per cent of homogeneity (Af = 1.01 and D2 = 2.97, red dotted linos). 
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Figure 9. Homogeneity scale Ru as a function of b 2 o-g(z) 2 , as predicted by the ACDM model, for different thresholds approaching 
homogeneity (10, 1 and 0.1 per cent from homogeneity, coloured curves from bottom to top), for A^(< r) (left) and D2(r) (right). The 
corresponding WiggleZ results are shown as error bars of corresponding colour (the errors are found using lognormal realisations). The 
6 2 <rg(z) 2 values of the data increase with redshift slice (so the data points from left to right arc from low to high rcdshift). Not all rcdshift 
slices have measurements at 0.1 per cent (blue) since the data do not reach this value in those slices. The definition of homogeneity used 
by previous authors is where the data comes within lcr of homogeneity - we show these scales as black diamonds. It is clear that this 
definition has much greater stochasticity than our definition, which fits a smooth curve to many data points. Indeed, this approach gives 
quite uninformative results in this plane and cannot be compared to the model prediction. 



ability distribution would be correlated, and we do not cor- 
rect for this. The D 2 {r) measurement is much less correlated 
(see Figs. [2] and [3]). Also, as we have explained, D 2 (r) is the 
most robust measurement of homogeneity, so we have chosen 
this for our likelihood analysis. 

The likelihood distribution on the data is simply given 
by the mean and variance of the lognormal realisations in 
each bin. These give the expected variance of a ACDM dis- 
tribution sampled with the WiggleZ selection function, so 
take into account both cosmic variance and shot noise. Here 
we assume they are Gaussian-distributed by virtue of the 
central limit theorem. However we also consider the true 
distributions provided by the lognormal realisations - this 
gives similar results but with larger errors, see Appendix El 

At each r, the data therefore provide a probability dis- 
tribution pr>[D2(r)], which we model as Gaussians, as: 



Pd[D 2 (t-) 



-(0 2 (r)-M) 2 /2<r 2 



2tt 



(26) 



where \i is the WiggleZ value of D2 at radius r and a is the 
root mean square variance given by the lognormal realisa- 
tions. See Fig. 1 101 for an illustration. 

We also expect there to be be a likelihood distribution 
on the D 2 (r) value we would measure for a perfectly homo- 
geneous distribution, due to cosmic variance and shot noise 
caused by our selection function. We can represent this by a 
likelihood distribution on the homogeneity scale - this would 
not be a simple delta-function at D 2 (r) = 3, but would have 
some spread. This distribution should be one-sided - that is, 
we don't expect to measure D 2 (r) > 3, only D 2 (r) 3, if we 
have a distribution that approaches homogeneity|j We might 
expect it to be represented by the variance in the random 



catalogues, which are essentially homogeneous distributions 
sampled with the WiggleZ selection function. However, the 
same sources of noise are also present in the lognormal re- 
alisations, so we would potentially double-count errors if we 
also used lognormal realisations to determine the likelihood 
distribution of the data. This means we cannot easily deter- 
mine the true variance in the value of D 2 (r) we would expect 
to measure for a homogeneous distribution, independently 
of the variance of the data. 

For this reason, we choose to find the likelihood distri- 
bution of the data reaching 1 per cent of homogeneity. That 
is, we assume the likelihood distribution on the homogeneity 
scale, pii[D 2 (r)], is a delta-function at D2 = 2.97: 



p H [D 2 (r)] = 5[D 2 (r)- 2.97]. 



(27) 



We can then construct the cumulative probability dis- 
tribution function P(Rh ^ r), which gives the probability 
that the homogeneity scale has been reached at or before 
scale r, from: 



P(Rh < r) 



Pd [Dz(r) 



£> 2 (r) 



PD[D 2 (r)]dD 2 (r) 



Ph{x)<1x dD2(r 



(28) 



That is, the probability of having reached homogeneity is 
the area of the likelihood distribution of the data that falls 
at or above D2 = 2.97. 

This cumulative probability is calculated at each scale 
r. We can then find the PDF for the homogeneity scale, 
p(Rh), from 



2 In some cases the WiggleZ data does fall below Af = 1 or above 
D2 = 3; this can be explained as the effect of shot noise introduced 



by the selection function rather than a physical effect, as shown 
in our comparison with the GiggleZ simulation in Section 16.21 
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Figure 11. Probability distributions for the scale of homogeneity, 
p(Rh), for WiggleZ galaxies in each of the four redshift slices. The 
homogeneity scale is defined as the scale where the data reaches 
values 1 per cent away from D2 = 3, i.e. D2 = 2.97. 



p(Rh) 



dP{R H < r) 
dr 



(29) 



The PDFs for the homogeneity scale, for WiggleZ galax- 
ies in each redshift slice, are shown in Fig. 1111 We have inter- 
polated between the data bins in order to obtain smoother 
PDFs. We find the most probable Rh values from the mean 
of the distributions. These are all between 70 and 81 h~ x 
Mpc, and are listed in Table [4] They represent the most 
probable scale at which the galaxy distribution reaches 1 
per cent of homogeneity. We also list the values found 
for the bias-corrected data, which give the most probable 
homogeneity scales for the matter distribution, assuming 
ACDM+WMAP. 



6 ROBUSTNESS OF HOMOGENEITY 
MEASUREMENT 

In this section we address several issues that could poten- 
tially influence our measurement of the homogeneity scale, 
and perform tests to ensure the robustness of our results. 
We base our tests on the 15-hr 0.5 < z < 0.7 region, which 
is the largest and most populated WiggleZ sub-region, but 
the results are applicable to the entire survey. 



6.1 Fractal model test of selection function and 
boundary effects 

A major potential source of bias in our results is the method 
used to correct for edge effects and the selection function of 
the survey. We have used 100 random catalogues to cor- 
rect each individual WiggleZ measurement, as described in 
Section [3.11 However, this method can potentially bias ho- 
mogeneity measurements, since it weights measurements by 
the volume of spheres of radius r included in the survey, and 
so a ssumes a homogeneous distribution outside the surve y 
(e.g. IColeman fe Pietronero|[l992l ; ISvlos Labini et al.ll2009r i. 
It is therefore important to check that this is not impos- 
ing any distortion in our measured correlation dimension, 
so producing a 'false relaxation' to homogeneity. 

To test this we apply our correction method to a 
range of fractal distributions of known correlation dimen- 
sio n. This has been done prev i ously by a number of w orks , 
e.g. iLemson fe Sanders! |l99ll ). [Provenzale et al.l |l994 ) and 
|Pan fe Coles! (|2002r ). This allows us to check that our method 
returns the correct input correlation dimension up to the 
largest scales we measure in the survey, and to quantify any 
distortion that may occur. 

We generate our fractal distributions using the j3- 
model, a simple self-similar cascading model (see e.g. 
ICastagnoli fe Provenzale]|l99ll ). This method starts with a 
cube of side Ln and splits it up into M smaller cubes of side 
Ln/n (we take n = 2, so M = 8). Each subcube is then 
assigned a probability p of surviving to the next iteration. 
This is repeated for a certain number of iterations k, and the 
resultant set of survived points is taken as the final distri- 
bution. In the limit of an infinite number of iterations, this 
produces a monofractal with a correlation dimension given 
by 



lim 



log(pM) fe logpAf 



lop 



log 



(30) 



Our procedure for using /3-models to test our analysis 
method is as follows: 

(i) We choose a range of D 2 values (2.7, 2.8, 2.9, 2.95 
and 2.97), and for each we generate 100 fractal galaxy dis- 
tributions, with a boxsize of (Ln /i~ 1 Gpc) 3 , where Ln is the 
length of the longest side of the WiggleZ 15-hr 0.5 < z < 0.7 
selection function grid (623.5 fe _1 Mpc). 

(ii) We then sample each distribution with the WiggleZ 
selection function for the 15-hr 0.5 < z < 0.7 region. We nor- 
malise the resulting distribution to give the same number of 
points as WiggleZ galaxies in this region. This gives us frac- 
tal mock catalogues, and we then measure D2 for these in the 
same way as for the WiggleZ data, correcting the counts-in- 
spheres measurements with random catalogues. This gives 
a result that is influenced by both the WiggleZ selection 
function, and our correction method. 
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Fig. 1121 shows the mean Af(< r) and D2(r) results for 
the different fractal distributions up to D2 = 2.95, with the 
WiggleZ selection function and correction method applied. 
Even up to D2 =2.95 they are clearly distinguishable from 
ACDM and the WiggleZ data. The D2 values measured are 
consistent with the input values up to at least 200 /i _1 Mpc, 
well above the homogeneity scales measured for the data. 
These results indicate that the WiggleZ selection function 
and correction method do not have a significant effect on the 
measured correlation dimension up to the scales we measure 
for homogeneity. 

It is noticeable that the size of the error bars gets 
smaller for models with larger Z?2- This is a real effect in the 
model; for larger D2, more boxes survive at each iteration, 
so there are a smaller number of possible configurations for 
the final distribution, resulting in lower variance for a box 
of a particular volume. 

To quantify how well we can exclude fractal models, we 
fit a line of constant D2 to each set of fractal data, over the 
range [80, SOOj/i^Mpc (shown in Fig.[[2]). This gives us the 
best-fit D2 value we would expect to measure for each fractal 
distribution over this range, taking into account bias from 
the selection function. We then find the formal probability of 
these values fitting the WiggleZ data. Doing this, we find we 
can exclude a fractal dimension of -D2M = [2.9,2.95,2.97] 
at the [19,6,4]-<t level. In other words, we can exclude fractal 
distributions with dimension D2(r) < 2.97 at over 99.99 per 

cent confidence on scales from 80 to 300 fc~ 1 Mpc. 

Ou r results agree w i th th ose o f iLemson fc Sanders! 

|l99ll ). iProvenzale et all j 19941 ) and |Pan fc Colesl (|2002l ) 
who also find that for samples on scales larger than the 
homogeneity scale, boundary corrections do not have a sig- 
nificant effect on the analysis. A further check would be to 
test different types of fractal model other than the /3-model, 
but we leave this for future work. We also note that we still 
assume an FRW metric in our fractal analysis; an improved 
consistency check would be to calculate the actual metric in 
these fractal models (see Section [7]), but this is beyond the 
scope of this paper. 



6.2 Comparison with the GiggleZ A-body 
simulation 

We have also compared our results with a ACDM cosmologi- 
cal A-body simulation. This allows us to check that our an- 
alytical ACDM+WMAP model (incorporating Kaiser and 
streaming models for redshift-space distortions) is consis- 
tent, over the relevant scales, with a full simulation including 
nonlinear effects. It also provides a further test of selection 
function effects, since we can show that homogeneity mea- 
surements of a ACDM distribution are not distorted when 
the WiggleZ selection function is applied. 

The GiggleZ (Giga-parsec WiggleZ) simulation (Poole 
et al. 2012, in preparation) is a suite of dark matter A-body 
simulations run at Swinburne University of Technology, de- 
signed for theoretical analyses of the WiggleZ dataset. It was 
run using a modified version of the A-body code Gadget-2 
<|Springel et al.ll200lT ) using a WMAP-5 cosmology. We use 
the main simulation, which has a volume of (1000ft _1 Mpc) 3 
and 2160 3 particles of mass 7.5 x W 9 h~ 1 M Q . 

Halo finding for Gi ggleZ was performed using Sub- 
find (|Springel et al.ll200ll ), which utilises a friends-of-friends 
(FoF) algorithm to identify coherent overdensities of parti- 



cles and a substructure analysis to determine bound over- 
densities within each FoF halo. For our analysis, we rank- 
order the resulting Subfind substructure catalogues by their 
maximum circular velocity (V max ) and select a contiguous 
subset of 250 000 halos (selected to yield a number density 
comparable to the survey) over a range of V max chosen to 
yield a bias comparable to that of WiggleZ. We use a cat- 
alogue at a redshift of z — 0.593, corresponding to the mid 
redshift of the WiggleZ 0.5 < z < 0.7 redshift slice. 

We add the effect of redshift-space distortions by shift- 
ing the positions of the haloes according to their line-of-sight 
peculiar velocities. That is, the comoving position of each 
halo relative to an observer, x, is shifted by a vector Ax, 



Ax 



x (1+f) 
[x| H(z) 



■«rad, 



(31) 



where w rac j = (x ■ v)/|x| is the radial velocity of the halo 
along the observer's line-of-sight, and we place the observer 
at the same coordinates relative to the GiggleZ box as for 
the WiggleZ selection function grid. 

We then calculate Af(< r) and -D2M using two different 
methods: 

(i) Using the full GiggleZ box. We correct the measure- 
ment using a random distribution within the same volume 
box, with 100 times the number of galaxies as the GiggleZ 
sample. 

(ii) Applying the WiggleZ 15-hr 0.5 < z < 0.7 selection 
function to GiggleZ. This creates a mock WiggleZ survey 
containing 10 830 galaxies. We correct the measurement us- 
ing the random catalogues used for the WiggleZ data. This 
allows us to see the effects induced purely by the selection 
function. 

The results are shown in Fig. 1131 The full GiggleZ 
dataset is very consistent with both the WiggleZ data and 
the ACDM model. The consistency with the model indi- 
cates that the implementation of redshift-space distortions 
in the model, described in Section 14.21 has a good level of 
accuracy, to scales as small as ~ 20/i _1 Mpc. The deviation 
of the smallest-scale bin is a resolution effect, since the Gig- 
gleZ catalogue is sparser than WiggleZ (with 10 830 galaxies, 
vs. 17 928 WiggleZ galaxies in the same region and redshift 
slice). The GiggleZ results both with, and without, the Wig- 
gleZ selection function are also consistent within the size of 
the WiggleZ error bars, showing that the selection function 
and correction method do not have a significant effect. It 
can be seen that the selection function causes a few data 
points in the M{< r) plot to go below 1, and in the D2(r) 
plot to go above 3. This shows that adding shot noise can 
produce this apparent unphysical effect, explaining why this 
is also seen in some of the WiggleZ results (Figs. OandEJ. 



6.3 Comparison of different correction methods 

To further demonstrate the robustness of our correction 
method, we illustrate the results obtained for two alternative 
correction methods - firstly, using only complete spheres for 
the counts-in-spheres measurements (the 'exclusion' method 
mentioned in Section but still correcting for the selec- 
tion function, and secondly, using no correction for the se- 
lection function. By using the exclusion method we do not 
have to deal with survey edge effects; however, since Wig- 
gleZ is not volume-limited, it is still necessary to use random 
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Figure 12. Fractal model comparisons of the measured correlation dimension from the WiggleZ 15-hr region, 0.5 < z < 0.7 redshift slice. 
The WiggleZ data (black error bars) and ACDM model (blue curve) are compared with several different /3-modcls with different fractal 
dimension (£>2 = 2.7,2.8,2.9 and 2.95), which have been sampled with the 15-hr selection function and analysed in the same way as 
the WiggleZ data (coloured error bars). The uncertainties shown are the error-in-the-mean of 100 fractal realisations. The input fractal 
dimensions are shown as dotted lines with corresponding colours. The best-fit D2(r) values fit over the range r = [80, 300] h' 1 Mpc arc 
shown as solid lines. 




Figure 13. Comparison of the GiggleZ A'-body simulation with WiggleZ, for the 15-hr region 0.5 < z < 0.7 redshift slice. The Af(< r) 
results arc shown on the left, and D2(r) on the right. The WiggleZ data is shown as black data points, and a ACDM model is shown 
in blue. The results for the full GiggleZ box are shown as the red crosses. The green crosses show the results for the GiggleZ simulation 
sampled with the WiggleZ 15-hr 0.5 < z < 0.7 selection function. The measured homogeneity scale Rjj is indicated for each. 



catalogues to correct for the selection function. That is, we 
recalculate Eq. [4] but with G equal to the number of Wig- 
gleZ galaxies at the centre of complete spheres of radius r. 
We show this result in Fig |14l for the 15-hr 0.5 < z < 0.7 re- 
gion, as red error bars, where the uncertainty is calculated 
using our lognormal realisations. We compare this to the 
result using our correction method (black error bars). The 
red error bars are consistent with the black error bars, but 
show more scatter and have higher noise, which increases 
for larger radius, reflecting the lower statistics where fewer 
spheres contribute to the measurement. The measurements 
must also be cut off at a lower radius, since not enough larger 
spheres fit inside the survey region. 

We can also compare the result without any selection 



function correction, which illustrates the necessity of cor- 
recting for having a non-volume-limited sample. However, 
WiggleZ contains holes in the angular coverage, which are 
independent of assumptions about completion, so we must 
still take these into account. We therefore normalise each 
N l (r) measurement by the volume within the selection func- 
tion included in that sphere. So we calculate: 

Gcomplete(r) . . . VSF.i X p 

complete spheres % 

where G comp ictc('") is the number of WiggleZ galaxies at 
the centre of complete spheres of radius r, Vsf,i is the vol- 
ume within the selection function of the ith galaxy, and 
p = nw /VsF,tot is the mean density of WiggleZ, i.e. the mim- 
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ber of WiggleZ galaxies divided by the total volume in the 
selection function, excluding holes. We show this in Fig |14l 
as purple error bars, where the uncertainty is again calcu- 
lated using our lognormal realisations (so the uncertainty 
assumes ACDM). Although there is a vertical offset in the 
M{< r) plot, caused by the selection function, the D2(r) re- 
sults are remarkably similar to those when correcting for the 
selection function, though again with more noise. There is 
still a clear transition towards values close to D2 — 3 on 
large scales. The offset in Af(< r) can be attributed to a se- 
lection effect: since we use only complete spheres we weight 
the measurement towards the central part of the redshift 
range, where the completeness is highest and so the number 
density of WiggleZ galaxies is highest. Since Dz(r) is the 
slope of JV(< r), it does not depend on the number density 
itself and so is more robust when summing over a varying 
selection function. 



7 DISCUSSION 

Our WiggleZ Af(< r) and D2(r) results show a very strong 
agreement with an FRW-based ACDM model. However, one 
of the strongest arguments against previous homogeneity 
measurements is that the method of correcting for survey 
selection functions, such as using random catalogues as we 
do, can distort the data, producing a 'false relaxation' to 
homogeneity. We have therefore tested this by applying the 
WiggleZ selection function to a range of inhomogeneous frac- 
tal models and a ACDM A^-body simulation. We have shown 
that there is no significant impact on our homogeneity mea- 
surement, up to at least 200 /i _1 Mpc. In addition, we have 
compared the results from our correction method to an anal- 
ysis using only complete spheres, with and without correct- 
ing for the selection function, and have shown that these 
are consistent, further demonstrating the robustness of our 
result. We can rule out fractals with fractal dimension up to 
D2 = 2.97, on scales between 80 and 300 h~ 1 Mpc, at over 
99.99 per cent confidence. 

We can also be confident that our result is robust to 
any assumptions in modelling the Wi ggleZ survey selec- 
tion function, since iBlake et all (|201(t ) showed that even 
extreme variations in modelling the angular completeness 
produce only ~ 0.5a shifts in estimates of the power spec- 
trum. Changes in the parameterisation of the redshift dis- 
tribution were shown to cause larger deviations, but only at 
scales > 200 ft -1 Mpc, which are well above the scales on 
which we measure homogeneity. 

Our result is a very good consistency check of ACDM. 
However, it is not independent of the assumption of the 
FRW metric. A complication for all homogeneity measure- 
ments is that we can only o bserve galaxies o n our light cone, 
and not on spatial surfaces. IMaartensl |201ll ') points out that 
it is therefore not possible to make a homogeneity measure- 
ment without making some assumptions, such as the FRW 
metric (to convert redshifts to distances) and the Cosmologi- 
cal Principle. Indeed, we must always assume an FRW-based 
model, ACDM in our case, to convert redshifts and angles 
to distance coordinates. This is also problematic consider- 
ing that inhomogeneities produce perturbations in the FRW 
metric, so can potentially distort distance measur ements by 
affecting the paths travelle d by light rays (e.g. IWiltshirg 
120091 ; iMeures fc Brunill2012i 'l. However, it seems highly un- 



likely that our measurements would so closely agree with 
FRW-based ACDM (in both the amplitude and shape of 
the A/"(< r) and -D2M curves), if the distribution were, ac- 
tually, inhomogeneous up to the largest scales probed, and 
we had incorrectly assumed an FRW metric. We have tested 
our method using a fractal model, and shown that this gives 
a completely different form of these curves. Our assumption 
of FRW also seems reasonable, since we know that distance 
measurements from Type la supernovae fit the Hubble dia- 
gram well to first order up to z ~ 1.4 jConlev et a.1.1 [20~lll : 
iDavis et al.ll201l[ ). so that any perturbations due to inho- 
mogeneities can only be a second- or higher-order effect. All 
this means that we can take our results to be a strong con- 
sistency check of FRW-based ACDM. 

It would be possible to test this further by mak- 
ing isotropy measurements in thin redshift shells, over a 
range of redshifts, without converting them to distances. 
We could also directly test the assumption of FRW, by 
calculating the non-FRW metrics of our fractal models 
and testing our analysis on these. We could also calcu- 
late the effects of metric perturbations for ACDM, includ- 
ing the effects of backreaction, as in e.g. iBucherti J2000h. 
iLi fc Schwarzl (|2007l ). IWiltshirel <|2007bh and Bchrcnd et al. 
12008). This would allow us to test the effect of incor- 
rectly assuming an FRW metric. Alternatively, several pos- 
sible consistency tests of homogeneity and the Copernican 
Princi ple that do not assume FRW have been suggeste d 
by e.g. IClarkson et all J2008l); IShafieloo fc Clarksonl <|2010r) ; 
IMaartensl l|201ll ); iHeavens et al. I (|201ll 1. However, we leave 
these suggestions for future work. 

An alternative way of defi ning the 'homoge neity 
scale' was r e cently suggested by iBagla et alj (|2008l ') and 
lYadav et alj (|2010f l. as where the measured fractal dimen- 
sion D q becomes consistent with the ambient dimension 
within the la statistical uncertainty, aAD q ■ They make pre- 
dictions for a ACDM model, by deriving an approximation 
for D q and aAD q , given a particular correlation function, 
and showing how these scale with sphere size. Using this, 
they predict that the 'true' homogeneity scale in ACDM is 
260 h" 1 Mpc. They also predict that D q and aAD q scale 
the same way with the correlation function, and so their 
definition of homogeneity does not change with bias or red- 
shift. Their definition is therefore beneficial, as it is robust 
to the tracer galaxy population. It is also not arbitrary, and 
indeed the scale above which the fractal dimension is consis- 
tent with homogeneity within cosmic variance is arguably a 
physically meaningful scale to define as homogeneous, since 
above this scale the distribution cannot be distinguished 
from a homogeneous one. However, as we have pointed out 
fSection l3.6|) . it is difficult to apply their definition to a real 
measurement, since their approximation for aAD q only ac- 
counts for the variance of the correlation function (in the 
limit of weak clustering) and shot noise, but ignores errors 
due to survey geometry and the selection function. These 
additional contributions mean that real errors will always 
be larger, and so will always measure a smaller homogeneity 
scale than the 'true' one. Since these errors will be different 
for different surveys, and cannot be separated out from the 
variance in the underlying correlation function, homogeneity 
measurements made in this way cannot be easily compared 
between different surveys, or with theory. The benefit of our 
method of defining the homogeneity scale, even if we have to 
make an arbitrary choice about the value of AD q we accept 
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Figure 14. Illustration of the result using different correction methods, for the 15-hr 0.5 < z < 0.7 region. Black error bars show the 
result using our correction method with random catalogues. The red error bars are obtained when only complete spheres are used in 
the counts-in-spheres measurement, while still correcting for the selection function. The purple error bars show the result for complete 
spheres, with no selection function correction. In each case the uncertainties are calculated using lognormal realisations. The red and 
purple error bars are shifted slightly to the right for clarity. 



for homogeneity, is that it can be used to easily compare the 
results from different surveys and with a theoretical model. 

A homogeneity scale below 100 /i _1 Mpc may also seem 
to contradict the fact that the correlation function has a 
known feature, the B aryon Acoustic Oscillation (BAO) peak 
at ~ 105 fe~ 1 Mpc (lEisenstein et all 12005 : IPercival et~ai1 
l20ld ; iBeutler et al.llioTll : iBlake et al"[|2011d ). However, the 
BAO peak has only a small impact on the counts-in-spheres 
statistic. It is visible in our ACDM prediction for D2(r) 
(which is more sensitive to it than J\f(< r), since it is a differ- 
ential measurement), as a small dip at just over 100 /i -1 Mpc 
(see Fig[7J). It means that the D2(r) curve does not increase 
monotonically around this scale, so we must be careful if 
attempting to measure an intercept with homogeneity that 
lies close to this. The magnitude of the distortion due to the 
BAO peak, AD2, is of order ~ 0.01 for b ~ 1 so does not af- 
fect our homogeneity scale measured at D2 = 2.97. However, 
it is more signi ficant for more hig hly-biased tracers (as also 
pointed out by iBagla et al.l 120081 ) , and for a highly-biased 
population, such as LRGs, it may be necessary to measure 
the homogeneity scale at a D2 value closer to 3, to avoid 
this region. 

It has been pointed out (e.g. ISvlos Labini fc Pietronerol 

I2010T I that measurements of large structures in galaxy sur- 
veys are seemingly at odds with a homogeneity scale be- 
low 100 ft. - Mpc. Previous survey s have detected struc- 
tures on scales much large r than this l|Geller k, Huchral 19891 ; 
Ide Lapparent et alj H9891. The largest observed structure in 
the Universe, the Sloan Great Wall, is 320 h~ x Mpc long 
(|Gott et al.ll2005h , appears inconsistent with the existence 
of a homogeneity scale below 100 ft -1 Mpc. However, this is 
not incompatible with our results, since we show the scale 
where the data has D2 ^ 2.97, and it is not impossible to 
have fluctuations on scales larger than this. Also, these large 
structures, including the Sloan Great Wall, are usually fila- 
mentary, whereas we have measured a volume statistic which 
averages over fluctuations. 



8 CONCLUSION 

We have measured the large-scale transition to homogeneity 
in the distribution of galaxies, using the WiggleZ survey, in 
four redshift bins between z = 0.1 and z = 0.9. We mea- 
sured the mean, scaled counts-in-spheres J\f(< r) and the 
correlation dimension, D2(r), and found these to be in ex- 
cellent agreement with a ACDM model with WMAP param- 
eters, including redshift space distortions. We also presented 
a new, model-independent method for determining the ho- 
mogeneity scale Rh from data. This involves fitting a poly- 
nomial curve to the data, and finding where this intercepts 
chosen values close to homogeneity. This is a more reliable 
method than finding where the data comes within la of ho- 
mogeneity, since it does not depend directly on the size of 
the error bars and is less susceptible to noise. It also allows 
a direct comparison between data and theory, and between 
different surveys of differing bias and redshift. 
We summarise our results as follows: 

• Our Af(< r) and D2 (r) results show a very strong agree- 
ment with a FRW-based ACDM+WMAP model incorpo- 
rating redshift-space distortions. They show a clear transi- 
tion from an inhomogeneous, clustered distribution on small 
scales, to a homogeneous one on large scales. This transi- 
tion matches that of the ACDM model. We have thereby 
conducted a very stringent consistency check of ACDM. 

• If we define the 'homogeneity scale' Rh as the scale 
where the data become consistent with homogeneity within 
1 per cent, then from a likelihood analysis of D2(r), we mea- 
sure R H to be 71 ± 8 /i _1 Mpc at z ~ 0.2, 70 ± 5 /i _1 Mpc 
at 2 ~ 0.4, 81 ± 5 h~ 1 Mpc at z ~ 0.6, and 75 ± 4 h _1 Mpc 
at z ~ 0.8. These values are consistent with those of the 
ACDM+WMAP model with best-fitting bias, of R H = 76, 
78, 81 and 78 7i _1 Mpc. 

• We find that the homogeneity scale of our 
ACDM+WMAP model increases with clustering am- 
plitude b(z)as(z). For a population with fixed bias, we 
therefore predict the homogeneity scale to grow over time, 
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since as(z) increases due to growth of structure. The bias 
of the WiggleZ galaxies increases with redshift, so the mea- 
sured Rh values do not change with redshift. If we correct 
our data for bias, assuming ACDM, then we measure a 
homogeneity scale for the matter distribution that increases 
over time, consistent with our ACDM+WMAP model. 

• The WiggleZ results are in excellent agreement with 
those of the GiggleZ A^-body simulation incorporating 
redshift-space distortions. It is also in excellent agreement 
with our analytic ACDM+WMAP model, showing that our 
model for redshift-space distortions is accurate down to 
scales as small as 20h~ Mpc. 

• We can exclude a fractal with fractal dimension up to 
D2 = 2.97, on scales between ~ 80/i _1 Mpc, and the largest 
scales probed by our measurement, ~ 300/t _1 Mpc, at 99.99 
per cent confidence. 

• By applying our analysis to the GiggleZ simulation, as 
well as a suite of fractal distributions of differing fractal di- 
mension, we have shown that our result is not significantly 
distorted by the WiggleZ selection function and our method 
of correcting for it. We also show that we obtain consistent 
results even using different correction methods, i.e. using 
only complete spheres for the measurement, with and with- 
out correcting for incompleteness. This therefore confirms 
the reliability and robustness of our results. 



9 ACKNOWLEDGEMENTS 

M.I.S. would like to thank Florian Beutler for many helpful 
comments on the paper. We also thank Eyal Kazin, David 
Parkinson, Luigi Guzzo and Andy Taylor for helpful dis- 
cussions and comments. M.I.S. acknowledges support from 
a Jean Rogerson Scholarship and a UWA Top-up Scholar- 
ship from the University of Western Australia. G.B.P. ac- 
knowledges support from two Australian Research Coun- 
cil (ARC) Discovery Projects (DP0772084 and DP1093738). 
The Centre for All-sky Astrophysics (CAASTRO) is an Aus- 
tralian Research Council Centre of Excellence, funded by 
grant CE11E0090. 

GALEX (the Galaxy Evolution Explorer) is a NASA 
Small Explorer, launched in April 2003. We gratefully ac- 
knowledge NASA's support for construction, operation and 
science analysis for the GALEX mission, developed in co- 
operation with the Centre National d'Etudes Spatiales of 
France and the Korean Ministry of Science and Technology. 

The WiggleZ survey would not have been possible with- 
out the dedicated work of the staff of the Australian Astro- 
nomical Observatory in the development and support of the 
AAOmega spectrograph, and the running of the AAT. 



REFERENCES 

Afshordi, N. 2004, Phys. Rev. D, 70, 083536, arXiv:astro- 
ph/0401166 

Amendola, L., & Palladino, E. 1999, ApJL, 514, LI, 

arXiv:astro-ph/9901420 
Bagla, J. S., Yadav, J., & Seshadri, T. R. 2008, MNRAS, 

390, 829, 0712.2905 
Baumann, D. 2009, ArXiv e-prints, 0907.5424 
Behrend, J., Brown, I. A., & Robbers, G. 2008, JCAP, 1, 

13, 0710.4964 



Beutler, F. et al. 2011, MNRAS, 416, 3017, 1106.3366 
Blake, C. et al. 2011a, MNRAS, 415, 2876, 1104.2948 

— . 2010, MNRAS, 406, 803, 1003.5721 

— . 2011b, MNRAS, 415, 2892, 1105.2862 

— . 2009, MNRAS, 395, 240, 0901.2587 

— . 2011c, MNRAS, 418, 1707, 1108.2635 
Blake, C, & Wall, J. 2002, Nat, 416, 150, arXiv:astro- 
ph/0203385 

Borgani, S. 1995, Phys. Rep., 251, 1, arXiv:astro- 
ph/9404054 

Brouzakis, N., Tetradis, N., & Tzavara, E. 2007, JCAP, 2, 

13, arXiv:astro-ph/0612179 
Buchert, T. 2000, General Relativity and Gravitation, 32, 

105, arXiv:gr-qc/9906015 
— . 2008, General Relativity and Gravitation, 40, 467, 

0707.2153 

Castagnoli, C, & Provenzale, A. 1991, A & A, 246, 634 

Clarkson, C, Bassett, B., & Lu, T. H.-C. 2008, Physical 
Review Letters, 101, 011301, 0712.3457 

Clifton, T., Ferreira, P. O, & Land, K. 2008, Physical Re- 
view Letters, 101, 131302, 0807.1443 

Coleman, P. H., & Pietronero, L. 1992, Phys. Rep., 213, 
311 

Coles, P. 1993, MNRAS, 262, 1065 

Coles, P., & Jones, B. 1991, MNRAS, 248, 1 

Conley, A. et al. 2011, ApJS, 192, 1, 1104.1443 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, 

AJ, 112, 839, arXiv:astro-ph/9606079 
Davis, T. M. et al. 2011, ApJ, 741, 67, 1012.2912 
de Lapparent, V., Geller, M. J., & Huchra, J. P. 1989, ApJ, 

343, 1 

Drinkwater, M. J. et al. 2010, MNRAS, 401, 1429, 
0911.4246 

Eisenstein, D. J. et al. 2005, ApJ, 633, 560, arXiv:astro- 
ph/0501171 

Ellis, G. F. R., & Buchert, T. 2005, Physics Letters A, 347, 

38, arXiv:gr-qc/0506106 
Fisher, K. B. 1995, ApJ, 448, 494, arXiv:astro-ph/9412081 
Fixsen, D. J., Cheng, E. S., Gales, J. M., Mather, J. C, 

Shafer, R. A., & Wright, E. L. 1996, ApJ, 473, 576, 

arXiv:astro-ph/9605054 
Gabrielli, A., Joyce, M., & Sylos Labini, F. 2002, Phys. 

Rev. D, 65, 083523, arXiv:astro-ph/0110451 
Geller, M. J., & Huchra, J. P. 1989, Science, 246, 897 
Glazebrook, K. et al. 2004, Nat, 430, 181, arXiviastro- 

ph/0401037 

Gott, III, J. R., Juric, M., Schlegel, D., Hoyle, F., Voge- 
ley, M., Tegmark, M., Bahcall, N., & Brinkmann, J. 2005, 
ApJ, 624, 463, arXiv:astro-ph/0310571 

Guzzo, L. 1997, NewA, 2, 517, arXiv:astro-ph/9711206 

Hamilton, A. J. S. 1992, ApJL, 385, L5 
— . 1993, ApJ, 417, 19 

Hatton, S., & Cole, S. 1998, MNRAS, 296, 10, arXiv:astro- 
ph/9707186 

Hawkins, E. et al. 2003, MNRAS, 346, 78, arXiv:astro- 
ph/0212375 

Heavens, A. F., Jimenez, R., & Maartens, R. 2011, JCAP, 

9, 35, 1107.5910 
Hogg, D. W., Eisenstein, D. J., Blanton, M. R., Bahcall, 

N. A., Brinkmann, J., Gunn, J. E., & Schneider, D. P. 

2005, ApJ, 624, 54, arXiv:astro-ph/0411197 
Joyce, M., Anderson, P. W., Montuori, M., Pietronero, L., 

& Sylos Labini, F. 2000, Europhysics Letters, 50, 416, 



© 2002 RAS, MNRAS 000, [TV?? 



20 M. I. Scrimgeour et al. 



arXiv:astro-ph/0002504 
Joyce, M., Montuori, M., & Sylos Labini, F. 1999, ApJL, 

514, L5, arXiv:astro-ph/9901290 
Kaiser, N. 1984, ApJL, 284, L9 

— . 1987, MNRAS, 227, 1 
Kolb, E. W., Matarrese, S., Notari, A., & Riotto, A. 2005, 

Phys. Rev. D, 71, 023524, arXiv:hep-ph/0409038 
Komatsu, E. et al. 2011, ApJS, 192, 18, 1001.4538 
Kurokawa, T., Morikawa, M., & Mouri, H. 2001, A & A, 

370, 358 

Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 
Lemson, G., & Sanders, R. H. 1991, MNRAS, 252, 319 
Lewis, A., Challinor, A., & Lasenby, A. 2000, ApJ, 538, 

473, arXiv:astro-ph/9911177 
Li, N., & Schwarz, D. J. 2007, Phys. Rev. D, 76, 083011, 

arXiv:gr-qc/0702043 
Maartens, R. 2011, Royal Society of London Philosophical 

Transactions Series A, 369, 5115, 1104.1300 
Martinez, V., & Saar, E. 2002, in Astronomical Data Anal- 
ysis II. Edited by Starck, Jean-Luc; Murtagh, Fionn D. 
Proceedings of the SPIE, Volume 4847, pp. 86-100 (2002)., 
ed. J.-L. Starck & F. D. Murtagh, Vol. 4847, 86-100, 
arXiv:astro-ph/0209208 
Martinez, V. J., & Coles, P. 1994, ApJ, 437, 550 
Martinez, V. J., Pons-Borderia, M.-J., Moyeed, R. A., & 
Graham, M. J. 1998, MNRAS, 298, 1212, arXiv:astro- 
ph/9804073 

Meures, N., & Bruni, M. 2012, MNRAS, 419, 1937, 
1107.4433 

Pan, J., & Coles, P. 2000, MNRAS, 318, L51, arXiv:astro- 
ph/0008240 

— . 2002, MNRAS, 330, 719, arXiv:astro-ph/0111234 
Peacock, J. A. 1999, Cosmological Physics, ed. Peacock, 
J. A. 

Peacock, J. A., & Dodds, S. J. 1994, MNRAS, 267, 1020, 
arXiv:astro-ph/9311057 

Peebles, P. J. E. 1980, The large-scale structure of the uni- 
verse 

— . 1993, Principles of Physical Cosmology (Princeton 

University Press) 
Percival, W. J. et al. 2010, MNRAS, 401, 2148, 0907.1660 
Pietronero, L., Montuori, M., & Sylos Labini, F. 1997, in 

Critical Dialogues in Cosmology, ed. N. Turok, 24-27, 

arXiv:astro-ph/9611197 
Provenzale, A., Guzzo, L., & Murante, G. 1994, MNRAS, 

266, 555 

Ras&nen, S. 2006, JCAP, 11, 3, arXiv:astro-ph/0607626 
— . 2011, Classical and Quantum Gravity, 28, 164008, 
1102.0408 

Ripley, B. D. 1977, Journal of the Royal Statistical Society. 

Series B (Methodological), 39, pp. 172 
Sachs, R. K., & Wolfe, A. M. 1967, ApJ, 147, 73 
Sarkar, P., Yadav, J., Pandey, B., & Bharadwaj, S. 2009, 

MNRAS, 399, L128, 0906.3431 
Scaramella, R. et al. 1998, A & A, 334, 404, arXiv:astro- 

ph/9803022 

Scharf, C. A., Jahoda, K., Treyer, M., Lahav, O., Boldt, E., 
& Piran, T. 2000, ApJ, 544, 49, arXiv:astro-ph/9908187 

Scherrer, R. J., & Weinberg, D. H. 1998, ApJ, 504, 607, 
arXiv:astro-ph/9712192 

Schwarz, D. J. 2002, ArXiv Astrophysics e-prints, 
arXiv:astro-ph/0209584 

Shafieloo, A., & Clarkson, C. 2010, Phys. Rev. D, 81, 



083537, 0911.4858 
Smith, R. E. et al. 2003, MNRAS, 341, 1311, arXiv:astro- 
ph/0207664 

Springel, V., Yoshida, N., & White, S. D. M. 2001, New A, 

6, 79, arXiv:astro-ph/0003162 
Sylos Labini, F. 2011, EPL (Europhysics Letters), 96, 

59001, 1110.4041 
Sylos Labini, F., Montuori, M., & Pietronero, L. 1998, 

Phys. Rep., 293, 61, arXiv:astro-ph/9711073 
Sylos Labini, F., & Pietronero, L. 2010, Journal of Statisti- 
cal Mechanics: Theory and Experiment, 11, 29, 1012.5624 
Sylos Labini, F., Vasilyev, N. L., & Baryshev, Y. V. 2009, 

A & A, 508, 17, 0909.0132 
Wiltshire, D. L. 2007a, New Journal of Physics, 9, 377, 

arXiv:gr-qc/0702082 
— . 2007b, Physical Review Letters, 99, 251101, 

0709.0732 

— . 2009, Phys. Rev. D, 80, 123512, 0909.0749 
Wu, K. K. S., Lahav, O., & Rees, M. J. 1999, Nat, 397, 

225, arXiv:astro-ph/9804062 
Yadav, J., Bharadwaj, S., Pandey, B., & Seshadri, T. R. 

2005, MNRAS, 364, 601, arXiv:astro-ph/0504315 
Yadav, J. K., Bagla, J. S., & Khandai, N. 2010, MNRAS, 

405, 2009, 1001.0617 
Zehavi, I. et al. 2005, ApJ, 630, 1, arXiv:astro-ph/0408569 



APPENDIX A: GAUSSIAN VS TRUE 
PROBABILITY DISTRIBUTIONS P D [D 2 (R)] 

In our likelihood analysis for the homogeneity scale (Sec- 
tion 15. 3p we assume that the distribution of the 100 lognor- 
mal realisations in each bin, po [-D2 (?")]; i s Gaussian. This 
allows us to interpolate between the data points and errors, 
so create more finely-spaced po[D2(r)] distributions, in or- 
der to determine a smoother PDF for the homogeneity scale, 
P(Rh ^ r). However, it is not obvious that these distribu- 
tions should be Gaussian. We therefore repeat the analysis, 
but use the true distributions given by the lognormal reali- 
sations. 

This gives the pd [D2 {r)] distributions shown in Fig. IA1I 
They are not smooth Gaussians, although they are close 
to Gaussian. Their resolution is limited by the number of 
lognormal realisations, so they could be improved by using 
more lognormal realisations, although we do not do this here. 

We then use these to calculate the PDF for the homo- 
geneity scale, P(Rh ^ r), in the same way as we did in Sec- 
tion [5]3] for the Gaussian distributions. This gives the PDFs 
shown in Fig lA2l for each redshift slice. The mean values 
and errors are shown in Table IA1I along with those from 
our original analysis. The values are very similar, though 
the errors are larger. This is because we cannot interpolate 
between data points as easily, and there are a finite number 
of lognormal realisations contributing to the distribution for 
each datapoint. This means the distribution is effectively 
smoothed, giving larger uncertainties. 
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Figure Al. The probability distributions PD[D2(r)] for each of 
the r bins (blue-to-red gradient indicates small to large radius) in 
the 0.5 < z < 0.7 redshift slice. Unlike in Fig llOl we do not assume 
these are Gaussian distributions; rather, we plot the distributions 
given by our 100 lognormal realisations. 
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Figure A2. Probability distributions for the homogeneity scale, 
p(Rjl), for WiggleZ galaxies in each of the four redshift slices. 
These are calculated from the probability distributions p£>[D2(r)] 
of the 100 lognormal realisations, rather than assuming Gaus- 
sians. 



Table Al. Comparison of the most probable Rh values from 
a likelihood analysis using the true po[D 2 (''')] distributions from 
lognormal realisations, and assuming Gaussian distributions. 



Redshift 


R H for true p D [D 2 (r)] 
distributions [h~ -"-Mpc] 


Rfj assuming Gaussians 
[fc-'Mpc] 


0.1 < 2 < 0.3 


79 ±19 


71 ±8 


0.3 < z < 0.5 


71 ± 13 


70 ±5 


0.5 < z < 0.7 


91 ± 16 


81 ±5 


0.7 < z < 0.9 


84 ± 16 


75 ±4 
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